Prepared for submission to JHEP 



Inflation with stable anisotropic hair: 
is it cosmologically viable? 



Sigbj0rn Hervik, a David F. Mota 6 and Mikjel Thorsrud 5 

"■Faculty of Science and Technology, University of Stavanger, N-4-036 Stavanger, Norway 
b Institute of Theoretical Astrophysics, University of Oslo, N-0315 Oslo, Norway 

E-mail: sigbjorn.hervik@uis.no, d.f.mota@astro.uio.no, 

mik j el . thorsrud@astro . uio . no 

Abstract: Recently an inflationary model with a vector field coupled to the inflaton was 
proposed and the phenomenology studied for the Bianchi type I spacetime. It was found 
that the model demonstrates a counter-example to the cosmic no-hair theorem since there 
exists a stable anisotropically inflationary fix-point. One of the great triumphs of infla- 
tion, however, is that it explains the observed flatness and isotropy of the universe today 
without requiring special initial conditions. Any acceptable model for inflation should thus 
explain these observations in a satisfactory way. To check whether the model meets this 
requirement, we introduce curvature to the background geometry and consider axisymmet- 
ric spacetimes of Bianchi type II, III and the Kantowski-Sachs metric. We show that the 
anisotropic Bianchi type I fix-point is an attr actor for the entire family of such spacetimes. 
The model is predictive in the sense that the universe gets close to this fix-point after a few 
e-folds for a wide range of initial conditions. If inflation lasts for N e-folds, the curvature 
at the end of inflation is typically of order ~ e~ 2N . The anisotropy in the expansion rate at 
the end of inflation, on the other hand, while being small on the one-percent level, is highly 
significant. We show that after the end of inflation there will be a period of isotropiza- 
tion lasting for ~ |iV e-folds. After that the shear scales as the curvature and becomes 
dominant around iV e-folds after the end of inflation. For plausible bounds on the reheat 
temperature the minimum number of e-folds during inflation, required for consistency with 
the isotropy of the supernova la data, lays in the interval (21,48). Thus the results ob- 
tained for our restricted class of spacetimes indicates that inflation with anisotropic hair 
is cosmologically viable. 
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1 Introduction 

It is plausible that inflation occurred around the energy scale of the grand unification. 
Since we do not have direct experimental access to such energy scales we need to be 
open-minded regarding the physics possibly occurring there. Indeed, when theoretical 
models of inflation are compared to observations they might provide the clue to understand 
fundamental physics at non-accessible energies. In this paper we are concerned with a 
specific model that violates isotropy; i.e. three-dimensional rotational invariance. Although 
rotational invariance is a well established feature of low-energy physics, several puzzling 
features of the large scale CMB anisotropics hints that this symmetry might have been 
broken during inflation [1-6]. At the same time recent progress in theoretical models has 
clarified that inflation with anisotropic hair is a theoretical possibility [7-16]. It has become 
clear, however, that specific realizations are often plagued by instabilities, either ghosts or 
unstable growth of the linearized perturbations [17-20]. 

In a series of recent papers M. Watanabe, S. Kanno and J. Soda studied the cosmology 
of a model with anisotropic hair which, apparently, is free of instabilities [21-25] (for 
related works by other authors also see [26-30]). Their model is inspired from supergravity 
which includes a massless vector field coupled to the scalar field (s). The vector part of 
the supergravity action has so far been neglected in cosmology, but for an approprately 
chosen coupling, the authors demonstrated interesting cosmological implications when the 
backreaction to geometry is properly accounted for. In particular, it is shown that the 
anisotropy in the expansion rate is stable and proportional to the slow roll parameter [21]. 
Even in the case when the anisotropy is very small, say on the micro level, the statistical 
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anisotropy imprinted in primordial fluctuations can be significant [24]. Thus, it has become 
clear that, in high precision cosmology, one cannot always neglect the back-reaction of 
vector-fields to geometry. 

One of the essential features of the conventional inflationary scenarios, however, is that 
they explain the observed flatness and isotropy of the universe today without requiring 
special initial conditions. For spatially homogeneous models of non-positive curvature 
containing comoving fluids obeying the strong energy condition, this is a consequence of the 
cosmic no-hair theorem [31-33] which guarantees that the curvature and shear of the spatial 
geometry decay rapidly during inflation. Given that inflation is stable for a sufficiently large 
time, the universe will, in agreement with observations, remain almost flat and isotropic 
until today. 1 The above mentioned papers, however, clearly demonstrate that the theorem 
may be violated if a massless vector field is appropriately coupled to the scalar field (and 
necessarily violating the assumptions given in the theorem). The phenomenology of the 
model has so far only been studied in the context of the spatially flat Bianchi type I model. 
The implications of curvature are therefore, until now, unexplored. A crucial question is 
whether the model, while violating the cosmic no-hair theorem, is still consistent with the 
observed isotropy and flatness of the universe today. It is therefore important to understand 
the behavior of the model with more general initial conditions. In this paper we shall 
introduce curvature to the model by considering axisymmetric spacetimes of Bianchi type 
II, III and the Kantowski-Sachs metric. To be specific we shall consider the well motivated 
case where both the potential of the scalar field and the coupling function between the 
vector and scalar fields are exponentials (this specific model is analyzed without spatial 
curvature in [23]). 

The organization of the paper is as follows. In section 2 we give a brief introduction 
to the model, while in section 3 we motivate our class of considered spacetimes. In section 
4 we derive the field equations. Section 5 is devoted to dynamical system analyses. First, 
in 5.1, we characterize the phase-space by identifying and classifying all fix-points. Then, 
in 5.2, we study the phase flow with arbitrary initial conditions. As we shall see the 
model is predictive in the sense that it provides unambiguous initial conditions for the 
post-inflationary era. Finally, in section 6, we show that these initial conditions provide a 
viable cosmological scenario. 

2 The model 

The starting point of the model is the following action for the metric g^, scalar field <f> 
and vector field A u : 



where (V0) 2 = g^V u 4> is the kinetic term of the scalar field, and = V ^Ay — VyA^ 
is the field strength of the vector field. In agreement with conventional notation g is the 

As a rule of thumb, 60 e-folds is sufficient. 
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determinant of the metric, R is the Ricci scalar and M p is the reduced Planck mass. 
We shall use metric signature (—1, 1, 1, 1). Motivated by dimensional reduction of higher 
dimensional theories we take the potential of the scalar field and the coupling functions to 
be exponentials: 

V(c/ ) ) = Voe X ^, (2.2) 
/(^) = / e Q 4, (2.3) 

where A and Q are constant parameters that characterize the model. We shall assume 
slow-roll inflation which implies A < 1. We shall treat the coupling constant Q as a 
free parameter, and study the implications on inflation for different values. We consider 
a positive potential which implies Vo > 0. There are no restrictions on the constant /o- 
Variation with respect to g^ v , (j) and A^, respectively, gives the equations of motion: 

Mffi = Ttf*) + , (2.4) 



V M T#*) = -Q^V,0, (2.5) 

V M T^) = +Q^V^, (2.6) 

where Et is the Einstein tensor, and Ca = —jf 2 (<j))F fll/ F IJ ' 1 ' is the Lagrangian for the vector 
field. Notice that the total energy-momentum tensor is conserved V^T^^ + T^j- A ^) = 0. 
It is manifest from the equations of motions that the coupling leads to exchange of energy 
and momentum between the scalar field and the vector field. The rate is determined by 
the coupling constant Q. It is the vector field that sources the shear degree of freedom of 
the expansion rate in this model. 2 Without the coupling to the scalar field, however, the 
vector field would decay rapidly leading to isotropization. 3 

Finally we write down the components of the energy-momentum tensors for the scalar 
field and the vector field: 

T#*> = Q(V0) 2 + V(<f>)\ + WV,0 (2.7) 
= -f*{<t>)F» a F« v - \f\4>)5>iF al 3F a ^ (2.8) 

where b„ =g m g av - 



3 On axisymmetric geometries 

The usual approach for anisotropic models where the shear is sourced by a homogeneous 
vector field, is to assume axisymmetric geometries, often called local rotational symmetric 

2 Since we do not consider any shear in the matter source, the term shear will unambiguously refer to 
the geometrical freedom in the expansion rate. 

3 As an interesting digression we mention that in a model where inflation is driven by non-Abelian gauge 
vector fields (without any scalar field), the anisotropic hair vanish exponentially fast although the field is 
sufficiently stable to drive inflation, see [34] and references therein. 
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(LRS) spacetimes. The vector field is assumed to be aligned parallel with the axis of 
symmetry, ie. orthogonal to the plane of rotational symmetry. Although this is a widely 
used assumption in the literature [10, 21-24, 35], it seems like a proper discussion is lacking. 4 
Beside being a simplifying assumption we think that the approach is appealing since, on 
the background level, one has the same symmetry in the spacetime geometry as in the 
(total) matter field. One can also loosely argue that if the symmetry was broken initially, 
spacetime would quickly become axisymmetric again, since there is no way to source a 
preferred direction in the plane orthogonal to the vector field. 5 The purpose of this section 
is to show rigorously that this is actually what happens in the Bianchi type I spacetime. 
We find that the universe isotropizes in the plane orthogonal to the vector field and that 
in the case of inflation this is a very rapid process. 

We consider the most general Bianchi type I geometry: 

ds 2 = -dt 2 + a 2 x dx 2 + a 2 y dy 2 + a 2 z dz 2 , (3.1) 

where a x (t), a y (t) and a z {t) are three independent scale factors. It is convenient to intro- 
duce three functions a(t), ai(t) and 02 (i) defined by 

a x = e a e- ai - a2 , a y = e a e a \ a z = e a e a2 . (3.2) 

In the special case cri(i) = 02 (i) we recover the axisymmetric geometry usually assumed. 

The three Hubble factors if, = — can then be written: 

1 a,i 

H x = a — &i — &2, H y = a + &i, H z = d + & 2 , (3.3) 

where a dot denotes differentiation with respect to the cosmic time t. It is also useful to 
introduce a mean Hubble rate 

H = ^{H X + H y + H z ) = a. (3.4) 

As in the Friedman-Lemaitre- Robertson- Walker (FLRW) metric, acceleration can be quan- 
tified in terms of the deceleration parameter q defined by: 

<? = -!- ^2- (3-5) 

Acceleration is then defined as q < 0. Note that q < is equivalent to ^e a > 0. The 
function e a can be interpreted as an isotropic scale factor. More precisely it is the geometric 
mean of the scale factors. Thus our definition of acceleration is equivalent to acceleration 



4 Ref. [24] comments that the shear variable (in the considered two dimensional hypersurface) is exponen- 
tially decaying in an expanding universe. But the physical interesting quantity is the dimensionless shear 
degree of freedom which, roughly speaking, measures the fraction between the anisotropic and isotropic 
parts of the expansion rate. As we shall see, this one is not necessarily exponentially decaying in an 
expanding universe. 

5 Note that in quadratic theories the curvature, in some sense, may source itself creating the possibility 
for non-axisymmetric Bianchi type I solutions [9]. 
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of the geometric mean of the scale factors. These definitions are standard in homogenous 
cosmologies, and generalize the definitions of the FLRW model in a natural way. 

To proceed we impose the metric (3.1) on the gravitational equation (2.4). We assume 
a single electric-type field. Since there are no spatial symmetries, apart from homogeneity 
of course, in the considered spacetime, we can without loss of generality align the field 
in the x-direction: F = ^F tll/ (t)dx lJ 'Adx u = Fio(t)dx Adt. In the gauge Aq = this field 
configuration corresponds to the vector potential A = A^dx^ 1 = A x {t)dx (where A x are 
related to F\q by F\q = —A x ). By taking linear combinations of the components of (2.4) 
one can eliminate the matter sources leading to the pure geometrical equation: 

(<t 2 — <7i) = — 3a(&2 — <ri )■ (3.6) 

We now introduce a shear degree of freedom: 

X± = E±zJh = (3.7) 
H a 

The quantity |-Xj_| is a measure of the anisotropy in the plane orthogonal to the vector 
field (as indicated by the subscript). For \X±\ = the hypersurface is isotropic and we 
have the axisymmetric metric usually assumed. Equation (3.6) then implies: 

^ = -\X ± \(2- q ). (3.8) 

In this equation we have changed to the number of e- folds, a, as time parameter by the 
identity ^ = H. From the gravitational equations it is easy to verify that q < 2 since we 
consider a positive potential V((f>) > 0. This implies -^g^ < 0. Thus we have shown that 
decays monotonically, although, in general, it will not necessarily decay rapidly. In 
the case of inflation, however, we have by definition q < 0. From (3.8) it then follows that 

d\ X I 

da < — 2|Xj_|. Thus inflation guarantees that |Aj_|(a) decays faster than the function 
e~ 2a . To summarize we have showed that the Bianchi type I spacetime isotropizes in the 
plane orthogonal to the vector field, and that in the case of inflation this is a very rapid 
(exponential) process. 

Let us finally use this result to prove that there will be no inflationary fix points with 
anisotropic expansion in the considered hypersurface. A dynamical system analysis for 
the spacetime (3.1) would require two shear variables, say X\ = anc l X2 = Note 
that from the definition of X± we have X± = X2 — X\. It follows that for a fix-point 
(^r = ^jjf = 0), we must have = 0. But, as shown above, this is not possible in 
our model. Thus there will be no additional fix-points if introducing a new shear degree of 
freedom. 

Although this analysis has been restricted to the Bianchi type I metric, we believe it 
sheds some light on the more general class of homogenous geometries, and we will use this 
as a motivation for the class of spacetimes to be introduced in the next section. 
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4 Field equation for a class of homogenous and axisymmetric spacetimes 



We shall now impose a class of homogenous spacetimes to the field equations. Motivated 
by the discussion above we consider a class of axisymmetric versions of Bianchi types I, 
II, III and the Kantowski-Sachs metric. We shall refer to these as BI, BII, Bill and KS, 
respectively. Type BI is spatially flat while BII, Bill and KS have anisotropic curvature 
in addition to the shear. For these spacetimes there is an intrinsic rotational symmetry in 
the spatial curvature which is aligned with the rotational symmetry of the expansion rate. 
Our considered class of spacetimes can be written on the form: 

ds 2 = _ dt 2 + e 2a(t) ( e -Mt) w lQ W l + e Mt) w ^ w 2 + e Mt) w ^ w ^ > (41) 

where w % are three time-independent and mutually orthogonal one-forms which are related 
to the coordinate basis in table 1. These spacetimes have rotational symmetry in the plane 
spanned by w 2 and w 3 . The rotational symmetry in the expansion rate is manifest from the 
line-element on this form. The rotational symmetry of the anisotropic curvature, on the 
other hand, is manifest first after an appropriate coordinate transformation. 6 To obtain 
the same symmetry in the (total) matter as in the spacetime geometry, we shall align the 
vector field parallel to w . 



Spacetime 


3 R 


w 1 


w 2 


w 3 


BI 





dx 


dy 


dz 


BII 




dx + k(ydz — zdy) 


dy 


dz 


Bill & KS 


2ke~ 2a ~ 2a 


dx 


{\-ky 2 Y x l 2 dy 


ydz 



Table 1. Geometric variables in the various spacetimes. w 1 are the one-forms in the line element 
(4.1). 3 R is the spatial Ricci scalar. A: is a constant. For Bill k < 0, while k > for KS. 



Since the one-forms are time-independent one can read the scale factors directly from 
the line-element (4.1): 

|| = e a - 2a , a ± = e a+CT , (4.2) 

where an and a± are the scale factors in the direction parallel and perpendicular to the 
vector field, respectively. The corresponding Hubble factors are: 

H\\ = a - 2a, H ± = a + a. (4.3) 

The mean Hubble rate becomes similar as in the previous section, H = a. 

We shall now introduce two spacetime dependent coefficients si and S2 with values 
specified in table 2. The former is defined by 3 R = —2(a + si&) 3 R, where 3 R is the three- 
dimensional Ricci scalar of constant time hypersurfaces. The functions 3 R for the various 
spacetimes are specified in table 1. The latter (52) determines the strength of a coupling 

6 As an example, in the Bianchi type II one can use y — rcos(^)) and z = rsin(0). 
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between the curvature and the energy density of the vector field. We can then treat the 
entire class of spacetimes in a unified way by expressing the equations in terms of si,S2 
and 3 R. 

Next, let us consider the matter sources. The most general field strength compatible 
with the class of spacetimes is on the form F = v w^/Kw 1 + (b + 2s2kv) w 2 Aw 3 , where 
v = v(t) is the dynamical degree of freedom and b is a constant. This field is homogeneous, 
axisymmetric and satisfy the identity dF = d 2 A = 0. Note that the field includes both 
an electric-type and a magnetic-type component which are parallel and pointing in the 
w 1 direction. Previous studies (which are restricted to the BI metric) have considered 
the case corresponding to b = 0, ie. neglecting the possibility of a magnetic component. 8 
Since our intention is to generalize the geometry, and not the matter fields, we will put 
6 = 0. In that case there is no magnetic-type component in the BI, Bill and KS spacetimes 
(since s 2 = 0). In the BII spacetime, however, where s 2 = 1, there is still a magnetic- type 
component. This is the minimal magnetic component required to satisfy the source- free 
Maxwell equations (dF = 0). Our considered field strength is therefore: 

F = v(t) w^w 1 + 2s 2 kv(t) w 2 Aw 3 . (4.4) 

The corresponding vector potential in the gauge Aq = is A = A^w^ = vw . 

We shall continue to specify all tensors, and perform all calculations, relative to the 
one- forms (dt,w l ,w 2 ,w 3 ). In this basis the energy- momentum tensor is diagonal and we 
write: 

T% = diag(-/>,p| |j pj.,p J .). (4.5) 

Since T% is a diagonal mixed tensor of rank (1,1), the components are invariant under a 
normalization of the one-forms (w l , w 2 , w 3 ), ie. under a change to an orthonormal basis. 
Thus the components of T% are physical quantities representing the energy density and 
pressure as measured in the fluid rest frame. We split the energy and pressure in the 
contributions from the scalar field (</>) and vector field (A): 

P = P4> + PA, 

P\\ =P<t> + (pa)||, (4.6) 
PL =P4> + (pa)±- 

The energy density and pressure of the scalar field takes the standard form in the entire 
family of spacetimes: 

P4> = l^ + V(<f>), P4> = \^-V(<P). (4.7) 

The energy density (and pressure) of the vector field, on the other hand, depends on the 
spacetime: 

PA = fe~ 2M Q^ 2 - s 2 ( 3 R)v 2 ^j . (4.8) 

7 As we shall see, it turns out that such a coupling exists in BII, but not in Bill or KS. 
8 Reference [22] treated magnetic fields perturbatively. 
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As S2 7^ only for BII it turns out that the spatial curvature couples to the energy density 
only in this case. Clearly, this is due to the magnetic-type field unique for BII. Note that, 
since 3 R < for BII, the coupling gives a positive contribution to the energy density. The 
equation of state, however, is similar for all spacetimes: 

(pa) || = ~PA, {pa)± = +PA- (4.9) 

The relation (jo^) || = — (p A )± also hold in the case of a massive vector field [28] (at least 
in the Bianchi type I metric). 

The field equations (2.4)-(2.6) for the considered family of spacetimes can be written: 



Spacetime 


Sl 


S2 


BI 








BII 


4 


1 


Bill & KS 


1 






Table 2. Spacetime dependent coefficients. 



1 1 2 3 7? 

3 R = -2(d + S!&fR, (4.13) 

p (f> + 3H(p (t> +p ( f > ) = QC A ^-, (4.14) 

26 

p A + 4(H + &)p A = -Q£ A ^, (4.15) 

where C A = \f 2 (4>) e ~ 2a+ia {v 2 + 2s 2 l v 2 ( i R)). From (4.12) it is manifest that the shear is 
sourced by anisotropic pressure and anisotropic curvature. 

We will now introduce dimensionless variables and rewrite the field equations as an 
autonomous set of first order differential equations. First we introduce the shear variable 

X=*^* = i, (4.16) 

and a variable for the curvature: 

«* - -sp- «■») 

Note that f2/< > in all spacetimes apart from KS where it is negative (and BI, of course, 
where it is zero). Furthermore, we need the variables: 

1 * e - a+2a v vfe- a+2a 
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The variables X, Y and Z are similar to those used in [23]. The additional variables £Ik 
and V are required to study the more general family of spacetimes considered here. We 
shall refer to the space spanned by the set of independent variables (X,Y, Z,Q,k,V) as 
state space. The space spanned by the constant parameters A and Q, we shall refer to as 
parameter space. Note that V^O only in BII due to the coefficient S2 in the definition of 
V. It is only in BII, where 3 R couples to the energy density of the vector field, that the 
variable V is needed. 

We can make use of the Hamiltonian constraint equation (4.10) to eliminate V{4>) from 
the equations of motions. Using the identity ^ = H we change to the scale, a, as time 
parameter, and the autonomous equations becomes: 



dX 
da 
dY 
da 



1 



Z 2 (X + 1) + X 



3(X 2 



1) + -Y 2 
1 2 



(Y + A) (s(x 2 - 1) + ^y 2 ) + \yz 2 + (Q+ ±)Z 



n K ( Sl + x + 4v 2 (i + x)) 

A, 



+ n K (y + 3A + 2V 2 (3A -6Q + 2Y)) , 



— =Z 
da 



3(X 2 



1) + -Y 2 

' 2 



da 
dV 
da 



20 



K 



QY + l 



1 



2X + -Z 1 
3 

1 



+ n K [Z + AZV 2 - 12V] 



1-81X + zx z + -y z + -z* + n K + asi k v 



(QY + 2X- 1)V + s 2 Z, 



The dynamical variables are subject to the constraint: 

x 2 + \y 2 + \z 2 + n K (1 + 2V 2 ) < 1, 

6 6 



(4.19) 
(4.20) 

(4.21) 
(4.22) 
(4.23) 

(4.24) 



which follows from our considered case of a positive potential (V((j)) > 0). For BI, BII 
and Bill this implies absolute upper bounds on each of the variables since £Ik > 0. More 
specifically X must be in the interval (—1, 1), Y and Z in \/6), while Qk (l + 2V 2 ) 

in (0,1). In the KS spacetime, however, each of the variables might be arbitrarily large 
individually since £1% < 0. Due to the negative curvature, the KS spacetime might collapse. 
At the turning point, from expansion to contraction, we have H= 0, and all the variables 
(apart from V) diverges (— > 00) as seen from the definitions (4.16)-(4.18). For initial 
conditions leading to collapse, a is then usually not a suitable time-parameter. For our 
purposes, however, it is fine since we are just interested in whether, and eventually when, 
the universe collapse for a given set of initial conditions. 

It is only in BII that both 3 R 7^ and s 2 7^ 0. For the other spacetimes the equations 
simplifies somewhat: 

n K -»• for BI, 

V -> for BI, Bill and KS. 



Although we shall study the dynamics in terms of the independent variables (X, Y, 
Z, Q-K, V), it is useful to introduce some auxiliary variables in order to make a closer 
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connection to the physics. The (energy) density parameters for the vector field and the 
scalar field can be expressed in terms of the independent variables in the following way: 

a ^3Mjm = l z2 + 2a " v2 < 425 » 

and 

^ = SMjlP =^~X 2 -\z 2 - n K (l + 2V 2 ). (4.26) 

Furthermore, we can split the latter one in the contributions from the kinetic and the 
potential energy: = ^l^+flv, where J\i n = \Y 2 . The Hamiltonian constraint equation 
(4.10) can then be written on the generic form: 

i = x 2 + n kin + n v + n A + n K . (4.27) 

It is also useful to express the deceleration parameter defined in (3.5) in terms of the 
independent variables: 

q = -1 + 3X 2 + \y 2 + \z 2 + n K (1 + 4V 2 ) . (4.28) 



From the definition of £Ik it follows that 

da 



2n K ( Sl X-q). (4.29) 



We notice that, due to the shear , there is no guaranty for monotonically decaying curvature 
during inflation (q< 0). Nevertheless, as we shall see in the next section, the model turns 
out to be predictive since the potential energy of the scalar field fly is monotonically 
increasing in a large region of state space. This leads the universe close to an anisotropic 
fix point with linear stability and vanishing curvature. 



5 Phase space analysis 

Equipped with the field equations we shall now investigate the phase-space structure by 
dynamical system analysis and simulations. First we shall identify the fix points of the 
system and classify their linear stabilities. Although this is a powerful way to characterize 
phase space qualitatively, it gives unambiguous predictions only in the linear regime close 
to the fix points. With arbitrary initial conditions it is therefore usually necessary to run 
simulations. As we shall see, however, in our case the potential energy of the scalar field is 
monotonically increasing in a large region of state space for BI, BII and Bill, leading the 
system close to a stable and unique anisotropic fix point of type BI. The KS spacetime is 
a bit more complicated. Firstly, in this case f2y is not a monotone function. Secondly, the 
KS spacetime might collapse. We therefore investigate the phase flow in the KS universe 
by running simulations. Our results indicates that if the universe does not collapse, its fate 
is similar to that of the Bianchi type spacetimes. The model is therefore predictive and, 
as we shall demonstrate by simulations, the universe gets close to the stable anisotropic 
fix-point typically within a few e-folds. 
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5.1 Fix-points and linear stability 

The fix points of the system are found by setting the left-hand side of the dynamical 
equations (4.19)-(4.23) equal to zero and solving the algebraic equation. The stability is 
determined by linearizing the field equations around the fix-points, ^p- = M5X 1 , and 
evaluating the eigenvalues of the matrix A4 . If the real part of all eigenvalues are negative, 
the fix-point is stable and we call it an attractor. If not all values are negative, the fix-point 
is unstable. Unstable fix-points are called saddles if there are both positive and negative 
eigenvalues, and repellers if all are positive. In general the phase flow goes from repellers, 
possibly via saddles, towards attractors. Without loss of generality we shall assume A > 0. 9 
Under the assumption A <C 1 we find 6 fix-points satisfying the constraint (4.24) or laying 
on its boundary. The fix-points are named (a)-(f) and the positions summarized in table 
3. We shall first give a brief overview before studying each of them more carefully. 

Name X Y Z 2 n K V 

~TT\ XQ-2 ^2 6AQ-12 n s 2 Z 

(b) -A 

(r \ _ q n 18A 3(2Q+A) n 

1 u 2Q-A 2Q-A u 

(r\\ 2-A 2 3A n 12-3A 4 n 

W _ 2+2l? ~T+A? U ~4(1+A 2 ) 2 U 

(e) free ±\/6 - 6X 2 

(f) I o f 

Table 3. Coordinates of fix-points. For fix-point (a) we have written down an approximation. The 
exact position of (a) can be found in the text. 



Name Spacetime Existence Stability 



(a) 

(b) 
(0 
(d) 
(e) 
(f) 



<1 



Comment 



BI 
FLRW 
KS 
KS 
BI 
Bill 



Q>1 
2Q > A 



attractor 

saddle 
saddle 
unstable 
saddle 



1 



■1 + 



3A 2 
2+2A 2 



2 
1/2 



boundary 
boundary 



Table 4. Properties of fix-points. The stability of (b) depends on A and Q (it is a saddle when (a) 
exists, and stable if not). See text for more details. 



In table 4 we give an overview of certain properties of the fix-points. Notice that the 
deceleration parameter q ~ — 1 for (a)-(d) while q > for (e) and (f). Thus (a)-(d) are 
inflationary fix-points, while (e) and (f) are decelerating. The fix-points (e) and (f) are on 
the boundary, ie. V(4>) = 0, while (a)-(d) satisfy the constraint (4.24). Notice that there are 

9 The situation A — > —A is equivalent to 4> — > —(f>. 
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no inflationary fix-points of type BII or Bill. For KS there are two inflationary fix-points, 
(c) and (d), but they are both saddles. A fix-point of special significance is the anisotropic 
attractor (a). Notice that it exists only when Q S> 1. Thus, interestingly, the shear X is 
small if it exists. The only fine tuning in the model is therefore the usual A <C 1 required 
for slow-roll inflation. No additional fine-tuning in any of the parameters is needed to avoid 
a dominating and observationally unacceptable shear. The model therefore represents a 
serious alternative to more conventional models. Fix-point (a) where identified in [23] and 
the stability determined within the Bianchi type I framework. Note, however, that (a) is 
on the border to the entire family of spacetimes as all becomes type BI in the limit of zero 
curvature. It is therefore important to check the stability with respect to a broader class of 
spacetimes. As we shall see, (a) turns out to be the unique attractor for the entire family 
of spacetimes. Loosely speaking, we can therefore say that the curved spacetimes becomes 
more and more Bianchi type I as they converges towards (a). Technically, spacetimes of 
type BII, Bill or KS never becomes BI, of course, but with decaying curvature they can 
come arbitrarily close. See figure 1 for a simulation close to the attractor (a). Note how all 
spacetimes converges towards a common point which represents fix-point (a). In figure 2 we 
show a simulation of type KS where the phase flow goes via the saddles (c) and (d) before 
converging towards (a). Fix-point (b) is isotropic and of type FLRW. In [23] it was shown 
that (b) is a saddle in the parameter region where (a) exists, and a stable attractor where 
(a) does not exist. As we shall see this holds also in the more general class of spacetimes 
considered here (apart from BII if Q <C — 1). 

We shall now examine each of the fix-points more carefully. 



Fixpoint (a) This is the anisotropically inflationary fix-point of Bianchi type I first 
identified in [23]. In table 3 the coordinates of (a) are given only to lowest order in the 
small quantities ^ and Q~ 2 . The exact position is: 



2(A 2 + 2QX - 4) 
~ A 2 + 8QA + 12Q 2 + 8' 
12(A + 2Q) 
A 2 + 8QA + 12Q 2 + 8' 
18(A 2 + 2QA - 4)(-A 2 + 4QA + 12Q 2 + 



z 2 = 



(A 2 + 8QA + 12Q 2 + 



£1k '■ 
V = 



0. 



s 2 Z 



1 - 2 X - QY' 



(5.1) 
(5.2) 

(5.3) 
(5.4) 
(5.5) 



Since the energy of the vector field is positive, ie. Qa > 0, we get the condition Z 2 > 0. The 
fix-point (a) therefore only exist in the parameter region where 2QX + A 2 > 4. Since A <C 1 
this implies Q> 1. As seen most directly from the approximations in table 3, it follows 
that (a) is near the origin of the state space (X,Y,Z,Qk,V)- Consequently the fix-point is 
strongly dominated by the potential energy of the scalar field. For clarity we expand the 
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eigenvalues in A and Q 2 , and truncate at zero order. For BI the three eigenvalues are: 



f -3, -| - iy 3(2AQ + A 2 - 4) - H -| + ^ 3(2AQ + A 2 - 4) - J I . (5.6) 

For BII one has the additional eigenvalues —3 and —2, while for Bill and KS one has the 
additional eigenvalue —2. The real part of all eigenvalues are negative for all spacetimes. 
Thus we have showed that the flat and anisotropic fix point identified in [23] , is an attractor 
also for the more general class of spacetimes considered here. Moreover, as we will verify 
below, it turns out to be the unique attractor (when it exists) for our considered class. Fi- 
nally we mention that, as shown in [23], there is an exact power-law solution corresponding 
to this fix-point, where the line element takes the form: 

ds 2 = _^ 2 + t 2k ^dx 2 + t 2fcl+2fc2 (dy 2 + dz 2 ), (5.7) 
where, to lowest order, k\ ~ ^ and k.2 ~ X %\Q . 

Fixpoint (b) This is a flat inflationary fix-point of type FLRW containing a scalar field 
dominated by its potential energy. The three eigenvalues for BI are: 

I(2AQ + A 2 -4), -3 + ^A 2 , -3 + ^A 5 

For BII one has the additional eigenvalues — 1 — XQ and —2 + A 2 , while for Bill and KS one 
has the additional eigenvalue — 2 + A 2 . This means that for the entire class of spacetimes, 
(b) is a saddle when (a) exists. If (a) does not exist, ie. 2XQ + A 2 — 4 < 0, then (b) is an 
attractor for BI, Bill and KS. For BII it is also stable in a large parameter region if (a) 
does not exist, but not if Q <C — 1, in which case it is unstable. 

Fixpoint (c) This is an inflationary fix-point of type KS containing an electric-type field 
and a cosmological constant (since the kinetic part of </> is vanishing). Essentially it is a 
generalization of the solutions with a pure cosmological constant found in [36] and [37]. 
From the condition Z 2 > it follows that the fix-point only exist in the parameter region 
2Q > X. The four eigenvalues are: 



Notice that the real part of the two latter eigenvalues are always negative when 2Q > A, 
ie. when it exists. We note that (c) is a saddle. 

Fixpoint (d) This is an inflationary fix-point of type KS containing only our considered 
scalar field. First note that there are no condition on Q for the existence of this fix-point. 
To lowest order in A the eigenvalues of (d) are: 

(-3 + 0(A 2 ), -6 + 0(X 2 ), 3 + 0(X 2 ), 3XQ + 0{X 2 
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Since A <C 1 the first and second eigenvalues are negative, while the third are positive. The 
sign of the last eigenvalue depends on Q. In any case (d) is a saddle. Although our analysis 
focus on slow-roll inflation A <C 1, we mention that there exist a region in parameter space 
(where A is larger than unity) where all the eigenvalues are negative and the fix point is 
a stable attractor. In this region, however, the deceleration parameter is positive and the 
fix-point is not inflationary. 

Fixpoint (e) This is a decelerating fix-point of type BI containing only a pure kinetic 
scalar field (thus the fluid is stiff p = p). The solution is part of a broader solution commonly 
referred to as Jacobs disc [38]. See also [39] for a discussion of such solutions. Note that 
(e) is a curve of fix-points, satisfying 6X 2 + Y 2 = 6. Since V(<p) = 0, on the curve, it lays 
on the boundary of our considered state space. We have two sets of eigenvalues depending 
on the sign of Y = ±y/6 — 6X 2 . The eigenvalues for BI: 

(o, 1-2XtQVQ-QX 2 , 6±A^6-6X 2 ) 

The sign of the third eigenvalue is always positive since A <C 1, while the second eigenvalue 
depends on Q and the sign of Y. For BII one has the additional eigenvalues 4 — 8X and 
-1 + 2X± QV6 - 6X 2 , while for Bill and KS one has the additional eigenvalue 4 — 2s±X. 
In any case (e) is unstable. 

Fixpoint (f ) This is a decelerating fix-point of type Bill containing no fluids. Essentially 
it is a Bianchi type III generalization of the Milne universe (but with a trivial flat direction). 
Like (e), it lays on the boundary V{4>) = 0. The four eigenvalues are: 

(3, -3/2, -3/2, -3/2). 

Thus (f) is a saddle. 

5.2 Phase flow with arbitrary initial conditions 

The stability analysis above determines the phase flow close to the fix-points. With ar- 
bitrary initial conditions, however, we need something more to determine the fate of the 
dynamical system. For the spacetimes of type BI, BII and Bill we shall now see that far 
away from (a) the phase flow can be characterized by the potential energy of scalar field. 

In section 4 we defined an auxiliary variable Oy, representing the potential energy of 
the scalar field. The equation of motion is: 

^ = Q v [XY + 2q + 2} = n v F, (5.8) 
da 

where F = GX 2 + (Y + ±) 2 -^- + lZ 2 + 2Q K (l+4V 2 ). Note that if Q K > and \X*\ > A for 
at least one of the variables, then F > 0, and consequently, £ly is monotonically increasing. 
Note that, as implied by (4.27), £ly ~ 1 when all the independent variables are small (<C 1). 
Regardless of initial conditions, the system will therefore approach some solution where 
Qy ~ 1 and all the independent variables are small (<S 1). The monotone function Qy 
will therefore lead the system close to the origin of the state space (X,Y,Z,Qk,V). In this 
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fl_K 0.002 




0.0010 

Figure 1. The phase flow of X, Z and SIk with A = 0.1 and Q = 50. The black, green, red and 
blue curves, respectively, represents simulations with BI, BII, Bill and KS initial conditions close 
to fixpoint (a). All trajectories converge towards a common point which is the type BI anisotropic 
fix point (a). 

region we have the stable anisotropic fix-point (a). The system will therefore eventually 
end up at the anisotropic fix-point (a) for arbitrary, but non-special 10 , initial conditions. 

In principle, it could take arbitrary long time for the system to get close to (a), since it 
might spend long time at the some of the saddles (for example at (b) which is also located 
close to the origin). Simulations, however, demonstrates that the typical time spent to get 
close to (a), starting with arbitrary initial conditions, only represents a minor fraction of 
the total number of e-folds during inflation. To demonstrate this, we define the distance 
to the fix-point (a) by 

d = y/(x - x a y + {y- Y a y + (z- z a y + (n K - (n K ) a y + (v - v a ) 2 , (5.9) 

where X a is the value of X at (a) and similar for the other coordinates. Note that the 
distance between (a) and (b) is of order ~ A ~ Q ■ We say that the system is close to (a) 
when D < ^. With this definition we ensure that the system is effectively unaffected by 
(b) and attracted by (a). In table 5 we have summarized some simulations with A = 0.1 
and two different values for Q. For Q = 50 we note that it typically takes 4-14 e-folds for 
the system to get close to fix-point (a). For Q = 1000 it typically takes 5-7 e-folds. 

10 One can construct special initial-conditions not leading to (a). With initially vanishing vector field for 
instance, Z = 0, the universe will converge towards the isotropic fix-point (b) instead of (a). 
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Figure 2. The phase flow of X, Y and Z for KS initial conditions with A = 0.1 and Q = 50. 
Initial conditions and direction of flow is indicated by the arrow. The fixpoint (a) , (b) , (c) and (d) 
are indicated by the colored points. Initial conditions are carefully chosen such that the solution 
"rides" on both KS-saddles (d) and (c) before ending up in the BI anisotropic attractor (a). 

The above analysis does not apply for KS since Sly is not a monotone function when 
£Ik < 0. Instead we study the KS spacetime by numerical simulations. Some results for 
randomly chosen initial conditions are summarized in table 6. As mentioned above, the KS 
universe might collapse, in which case H — > 0. Simulations shows that non-special initial 
conditions, either lead to collapse or the system will converge towards the fix-point (a). For 
initial conditions leading to collapse, the table shows how long time it takes before H = 0. 
We note that H — > very quickly, typically within the first third of the first e-fold, if the 
universe collapses. For initial conditions not leading to collapse, simulations shows that 
the system ends up at (a), and the table shows how long time it takes to get a distance 
D < ^ from (a). The results are very similar as those for BI, BII and Bill if the universe 
does not collapse. The system is close to (a) typically after 4-14 e-folds for Q = 50 and 5-7 
e-folds for Q = 1000. 

Supported by simulations we have showed that for initial conditions leading to in- 
teresting cosmologies (no immediate recollapse), the universe gets close to (a) relatively 
quickly. 11 When the system is close to (a) we see from (4.22) that the curvature decays ex- 

11 We should mention that for several of the (randomly chosen) initial conditions in tables 5 and 6, the 
universe is not accelerating, ie. q might be positive at the initial time a = 0. The time taken to get close 
to (a) after the start of inflation, is therefore somewhat smaller than the times in the table for the cases 
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ponentially, tlx ~ e 2q . If inflation lasted for N e- folds, a rough estimate for the curvature 
at the end of inflation is therefore (rJx)end ~ e~ 2N . 

6 Late time behavior 

To check if the model gives rise to an acceptable cosmology we will now investigate the late- 
time behavior. We showed above that at the end of inflation the universe will be extremely 
close to the anisotropic fix-point (a), with a curvature of order ($7^- ) en d ~ e~ 2N if inflation 
lasted for N e- folds. At the end of inflation we assume a period of reheating where the 
energy of the scalar and vector fields are dumped primarily into radiation. Reheating is a 
poorly understood process, but we shall assume that it is sufficiently fast that we can take 
the initial conditions for the shear (X) and the curvature (£Ik) t° be similar as at the end 
of inflation. In our analysis we shall, for simplicity, introduce a single perfect fluid with 
equation of state p = up. The radiation dominated era after inflation can then be studied 
by setting u = ^, while the subsequent period dominated by non-relativistic fluids can be 
studied by setting u = 0. We shall neglect the cosmological constant, since it would not 
change the estimates we seek significantly. 

The conservation equation in the class of spacetimes yields 

p + 3H(1 + u)p = 0. (6.1) 

The associated density parameter is defined 

The Hamiltonian constraint equation can then be written: 

i = x 2 + n M + n K , (6.3) 

where X and is defined in the same way as above. Due to the constraint (6.3) the 
system is effectively two dimensional, and we choose Qm as the auxiliary variable. The 
autonomous equations for the post-inflationary era with a single perfect fluid can then be 
written: 

dX 1 

— = --X [(1 + 3oj)fl K + 3(1 - u)(l - X 2 )] + Sl Q K , (6.4) 
= Q K [(1 - n K )(l + 3w) - 2siX + 3(1 - uj)X 2 ] . (6.5) 

da 

This system has four fix-points (P1)-(P4) summarized in table 7. The real parts of the 
eigenvalues can be found in table 8. Note that there are no fix-points of type KS. Without 
introducing a dark energy in the model, the KS spacetime will eventually collapse [40]. 
The Bianchi type I fix-point (P4) represents the two special points on the Kasner circle 
which are LRS [41]. Also note that it is a special case of (e) in table 3. Setting lo = 0, the 



where q > initially. In any case, since we typically assume inflation lasted for at least 60 e-folds, only 
minor fraction will be spend far away from (a). 
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attractors (P2) and (P3) represents the late time solutions for BII and Bill, respectively, 
for a universe without dark energy. Fix-point (P2) corresponds to the Collins-Stewart 
Bianchi type II exact solution [42], while (P3) is similar as fix-point (f) in table 3. For 
u = 0, (P2) and (P3) are the global attractors for the most general perfect fluids of type 
II and III, respectively, even including tilt [43], [44]. 

Note that the shear in the attractors (P2) and (P3), with X around unity, is quite 
extreme. Supernova la data gives the bound —0.012 < X to da y < 0.012 to one sigma 
confidence level for a Bianchi type I spacetime with rotational symmetry [45]. The goal of 
this section is to find the minimum number of e- folds during inflation, -/V m i n , required for 
a shear today in agreement with supernova observations. For this estimate we shall use 
l^todayl < 0.01. More stringent bounds can probably be found from the CMB which is 
very sensitive to shear [46-49]. For simplicity we shall use the supernovae since the final 
result is not very sensitive to the bound on X tQ( ^ ay . Thus, observational bounds certainly 
implies that the universe must have been very close to the flat FLRW saddle (PI) just 
before dark energy became significant. To determine -/V m i n we must therefore investigate 
the dynamics close to the fix-point (PI). Linearizing (6.4)-(6.5) around (PI) yields: 

£_ = -±(l- u )X + ai Q K , (6.6) 
da 2 

^ = (l + 3oj)n K . (6.7) 
da 

These equations can be solved exactly: 



3 



X(a) = e~s 



(1— uj)a 



(6.8) 



"5 + 3a; v ' 3Q 2 

fiA-(a) = ± e (1+3a;)Q " 27V , (6.9) 

where we have used the initial conditions X(0) = A ^, 2 and Qk(0) = ±e~ 2N . Here N is 
the number of e-folds during inflation. The initial condition for the shear corresponds to 
the shear at the anisotropic fix-point (a) , while for the curvature we have used the estimate 
at the end of inflation derived above. To get some intuition for the solution we shall now 
set u = 1/3 and study the behavior at different times. In the regime < a < §iV we 
have the approximation X ~ x ^ 2 e~ a . Thus there is an era of isotropization starting 
right after reheating, where the shear decays as ~ e _Q , lasting until a ~ |./V. Also notice 
that X ~ Qk around a = |iV. This is also true for KS since the sign of X will change 
from positive to negative around the same time. In the period |iV < a < N we have the 
approximation X ~ ±^sie 2a ~ 2N . In this period we have hi|X| ~ ln|f^| ~ 2a — 2N. 
We therefore say that the curvature "tracks" the shear in this period. Around a = N 
the parameters X and £Ik will approach unity and the universe is no longer close to the 
matter dominated fix-point (PI). The solutions (6.8)-(6.9) are then no longer valid. In 
the regime a > N the universe will in general converge towards the attractors (P2) and 
(P3), respectively, for BII and Bill, while KS will eventually collapse. Simulations of the 
autonomous equations (6.4)-(6.5) verifies these approximations. For an example with Bill 
spacetime, see figure (3). Note how the shear decays until it catches up with the curvature 
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Figure 3. Simulation of Bill after inflation with equation of state u> = 1/3 for the perfect fluid. 
The initial conditions comes from an inflation model with A = 0.1 Q = 50 and number of e-folds 
N = 60. 



slightly before a = |iV. After that the shear is tracked by the curvature, approaching unity 
around a ~ N. This behavior is in close agreement with the above approximations. It is 
easily seen from the solutions (6.8)-(6.9) that the same tracking behavior occurs also with 
uj 7^ 1/3. Thus the shear will be tracked by the curvature also in the matter dominated 
era. 

Reheating is complete when practically all of the energy is in radiation at thermal 
equilibrium. We shall denote the time and temperature at this stage by a rc h and T re h, 
the reheat time and temperature, respectively. Furthermore we shall denote the time at 
radiation-matter equivalence by a eq , and today by ao- Setting the initial time a rc h = at 
T = T rc h we get ao = ln-^ 3 -, where To = 2.73K is the temperature of the CMB today. 
Here we have used that the temperature of a photon gas scales as oc e~ Q . Matter-radiation 
equivalence occurred around redshift 3300 corresponding to a cq ~ ao — 8. We can then 
find the minimum number of e-folds during inflation, N m i n , required for a shear today 
consistent with supernova data, ie. |Ao| < 10~ 2 . The initial condition for the curvature is 
Qx(«rch) — e~ 2N . Assuming that the universe is close to the saddle (PI) all the way up to 
today, the curvature grows approximately as Q,k oc e 2a in the radiation dominated era, and 
as Qk oc e a in the matter dominated era. The curvature today is then (£Ik)o ~ e a cq+ a o-2N _ 
Since the shear is tracked by the curvature both in the radiation and matter dominated 
eras, we have Xq ~ (£Ik)o- Let |X max | denote the maximum shear today consistent with 
observations. The condition \Xq\ < \X max \, then leads to: 

iVmin = In ^ - I In \X max \ - 4. (6.10) 

As discussed above we shall use |X max | = 10 -2 , consistent with the supernova la data. A 
lower bound on the reheat temperature T re h > lMeV comes from the requirement that 
reheating must have occurred before the start of nucleosynthesis. This corresponds to 
-Wmin = 21. There is also an upper limit T^-g^ < 10 9 GeV if supersymmetry exists [50], 
corresponding to N^-n = 48. 
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To summarize, we have showed that plausible bounds on T Te ^ implies -/V m i n in the 
interval (21,48). 



7 Conclusion 

In this paper we have studied an inflationary scenario with a stable anisotropic hair. Since 
the model provides a counter-example to the cosmic no-hair theorem, it is important to 
study the model with more general initial conditions. Ultimately, as suggested in [21, 23], 
this might lead to a modified cosmic no-hair theorem. As a first step we took a dynamical 
system approach and showed that the stable anistropic fix-point of type BI, identified in 
[23], is the unique attractor for a wider class of spacetimes exhibiting spatial curvature. 
Moreover, for the considered Bianchi type spacetimes, we showed that the potential energy 
of the scalar field (fiy) is monotonically growing in a large variable region leading the 
universe close to the origin of the state space. In this region we have the stable anisotropic 
fix-point. This observation explains the rapid convergence in numerical simulations with 
arbitrary initial conditions, and, in fact, it is a reminiscent of the cosmic no-hair theorem 
itself. 12 The KS spacetime, however, which (depending on initial conditions) might col- 
lapse, is more complicated and our analysis relies on numerical simulations. For initial 
conditions not leading to collapse, our simulations shows that the behavior is quite simi- 
lar as for the considered Bianchi type spacetimes. For arbitrary, but non-special 13 , initial 
conditions, the universe will typically get close to the attractor after a few e- folds for any 
of the considered spacetimes. Thus, for the major fraction of e-folds required for inflation, 
the inflationary universe is well described by this fix-point. 

Thanks to these features the model provides unambiguous initial conditions for the era 
after inflation. After reheating, when the energy of the scalar and vector fields are dumped 
primarily into radiation, the only hair left is the anisotropic expansion. Interestingly, these 
initial conditions yields an acceptable late-time cosmology. After reheating there will be a 
period of isotropization lasting for ~ |iV e-folds, where N is the number of e-folds during 
inflation. After that the shear scales as the curvature and becomes dominant around N 
e-folds after the end of inflation. For plausible bounds on the reheat temperature, the 
minimum number of e-folds during inflation, required for consistency with the isotropy of 
the supernova la data, lays in the interval (21, 48). As a rule of thumb, successful inflation 
models must be stable for at least 60 e-folds. It is clear that this number is sufficient also 
for the considered model despite the stable anisotropic hair. So, to answer the question 
in the title: yes, the results obtained for our restricted class of spacetimes indicates that 
inflation with anisotropic hair is cosmologically viable. 



The only difference is that, while the no-hair theorem relies only on the fact that flv is a monotone 
function, Qv is only monotone until it is almost 1 in our considered model. Our argument therefore also 
relies on the dynamical system analysis. 

13 For example, if the energy density of the vector field is initially vanishing, it will always be vanishing, 
and the universe goes to the isotropic FLRW saddle point. 



- 20 - 



Acknowledgments 



DFM thanks the Research Council of Norway FRINAT grant 197251/V30. DFM is also 
partially supported by project PTDC/FIS/111725/2009 and CERN/FP/116398/2010. We 
also appreciate useful comments from an anonymous referee. 

References 

[1] G. Efstathiou, A Maximum Likelihood Analysis of the Low CMB MulUpoles from WMAP, MNRAS 
348 (2004) 885 [arXiv:astro-ph/0310207]. 

[2] A. de Oliveira-Costa, M. Tegmark, M. Zaldarriaga and A. Hamilton, The significance of the largest 
scale CMB fluctuations in WMAP, Phys. Rev. D 69 (2004) 063516 [arXiv:astro-ph/0307282]. 

[3] C. J. Copi, D. Huterer and G. D. Starkman, Multipole Vectors-a new representation of the CMB sky 
and evidence for statistical anisotropy or non-Gaussianity at 2 <= I <= 8, Phys. Rev. D 70 (2004) 
043515 [arXiv:astro-ph/0310511]. 

[4] H. K. Eriksen, F. K. Hansen, A. J. Banday, K. M. Gorski and P. B. Lilje, Asymmetries in the CMB 
anisotropy field, Astrophys. J. 605 (2004) 14, [Erratum-ibid. 609 (2004) 1198] 
[arXiv:astro-ph/0307507]. 

[5] F. K. Hansen, A. J. Banday and K. M. Gorski, Testing the cosmological principle of isotropy: local 
power spectrum estimates of the WMAP data, MNRAS 354 (2004) 641 [arXiv:astro-ph/0404206]. 

[6] N. E. Groeneboom, M. Axelsson, D. F. Mota and T. Koivisto, Imprints of a hemispherical power 
asymmetry in the seven-year WMAP data due to non-commutativity of space-time, 
arXiv:astro-ph / 101 1 .5353 . 

[7] L. H. Ford, Inflation Driven By A Vector Field, Phys. Rev. D 40 (1989) 967. 

[8] J. D. Barrow and S. Hervik, Anisotropically Inflating Universes, Phys. Rev. D 73 (2006) 023007 
[arXiv:gr-qc/0511127]. 

[9] J. D. Barrow and S. Hervik, On the evolution of universes in quadratic theories of gravity, Phys. 
Rev. D 74 (2006) 124017 [arXiv:gr-qc/0610013]. 

[10] L. Ackerman, S. M. Carroll and M. B. Wise, Imprints of a Primordial Preferred Direction on the 
Microwave Background, Phys.Rev.D 75 (2007) 083502 [arXiv:astro-ph/0701357]. 

[11] A. Golovnev, V. Mukhanov and V. Vanchurin, Vector Inflation, JCAP 0806 (2008) 009 
[arXiv:astro-ph/0802.2068]. 

[12] S. Kanno, M. Kimura, J. Soda and S. Yokoyama, Anisotropic Inflation from Vector Impurity, JCAP 
0808 (2008) 034 [arXiv:hep-ph/0806.2422]. 

[13] T. Koivisto and D. F. Mota, Vector Field Models of Inflation and Dark Energy, JCAP 0808 (2008) 
021 [arXiv:astro-ph/0805.4229]. 

[14] C. Germani and A. Kehagias, P-nflation: generating cosmic Inflation with p-forms, JCAP 0903 
(2009) 028 [arXiv:astro-ph/0902.3667]. 

[15] J. D. Barrow and S. Hervik, Simple Types of Anisotropic Inflation, Phys. Rev. D 81 (2010) 023513 
[arXiv:gr-qc/0911.3805]. 

[16] T. S. Koivisto and D. F. Mota, CMB statistics in noncommutative inflation, JHEP 1102 (2011) 061 
[arXiv: astro-ph /1011.2126]. 

[17] B. Himmetoglu, C. R. Contaldi and M. Peloso, Ghost instabilities of cosmological models with vector 
fields nonminimally coupled to the curvature, Phys. Rev.D 80 (2009) 123530 
[arXiv:astro-ph/0909.3524]. 



- 21 - 



[18] B. Himmetoglu, C. R. Contaldi and M. Peloso, Instability of anisotropic cosmological solutions 
supported by vector fields, Phys.Rev.Lett. 102 (2009) 111301 [arXiv:astro-ph/0809.2779]. 

[19] B. Himmetoglu, C. R. Contaldi and M. Peloso, Instability of the ACW model, and problems with 
massive vectors during inflation, Phys. Rev.D 79 (2009) 063517 [arXiv:astro-ph/0812.1231]. 

[20] S. M. Carroll, T. R. Dulaney, M. I. Gresham and H. Tarn, Instabilities in the Aether, Phys.Rev.D 79 
(2009) 065011 [arXiv:hep-th/0812.1049] 

[21] M. Watanabe, S. Kanno and J. Soda, Inflationary Universe with Anisotropic Hair, Phys. Rev. Lett. 
102 (2009) 191302 [arXiv:hep-th/0902.2833]. 

[22] S. Kanno, J. Soda and M. Watanabe, Cosmological Magnetic Fields from Inflation and 
Backreaction, JCAP 0912 (2009) 009 [arXiv:astro-ph/0908.3509]. 

[23] S. Kanno, J. Soda and M. Watanabe, Anisotropic Power-law Inflation, JCAP 1012 (2010) 024 
[arXiv:hep-th/1010.5307]. 

[24] M. Watanabe, S. Kanno and J. Soda, The Nature of Primordial Fluctuations from Anisotropic 
Inflation, Prog. Theor. Phys. 123 (2010) 1041 [arXiv:astro-ph/1003.0056]. 

[25] M. Watanabe, S. Kanno and J. Soda, Imprints of anisotropic inflation on the cosmic microwave 
background, MNRAS: Letters, 412 (2011) Issue 1, pp. L83-L87 [arXiv:astro-ph/1011.3604]. 

[26] R. Emami, H. Firouzjahi, S. M. Sadegh Movahed and M. Zarei, Anisotropic Inflation from Charged 
Scalar Fields, JCAP 1102, 005 (2011) [arXiv:astro-ph/1010.5495]. 

[27] R. Emami and H. Firouzjahi, Issues on Generating Primordial Anisotropics at the End of Inflation, 
arXiv:astro-ph/1111.1919. 

[28] K. Dimopoulos and J. M. Wagstaff, Particle Production of Vector Fields: Scale Invariance is 
Attractive, Phys. Rev. D 83 (2011) 023523 [arXiv:hep-ph/1011.2517]. 

[29] K. Dimopoulos, G. Lazarides and J. M. Wagstaff, Eliminating the n-problem in SUGRA Hybrid 
Inflation with Vector Backreaction, arXiv:astro-ph/1111.1929. 

[30] T. Q. Do, W. F. Kao and I. C. Lin, Anisotropic power-law inflation for a two scalar fields model, 
Phys. Rev. D 83 (2011) 123002. 

[31] R. M. Wald, Asymptotic behavior of homogeneous cosmological models in the presence of a positive 
cosmological constant, Phys.Rev.D 28 (1983) 2118. 

[32] I. Moss and V. Sahni, Anisotropy In The Chaotic Inflationary Universe, 
Phys. Lett. B 178 (1986) 159. 

[33] Y. Kitada and K. Maeda Cosmic no-hair theorem in power-law inflation, 
Phys.Rev.D 45 (1992) 1416. 

[34] A. Maleknejad, M.M. Sheikh-Jabbari and J. Soda, Gauge-flation and Cosmic No-Hair Conjecture, 
arXiv: hep-th/1109.5573v2. 

[35] C. G. Boehmer and D. F. Mota, CMB Anisotropics and Inflation from Non-Standard Spmors, Phys. 
Lett. B 663 (2008) 168 [arXiv:astro-ph/0710.2003]. 

[36] A.D. Linde and M.I. Zelnikov, Inflationary universe with fluctuating dimension, 
Phys. Lett. B 215 (1988) 59. 

[37] J. Yearsley and J.D. Barrow, Cosmological Models of Dimensional Segregation, 
Class. Quantum Grav. 13 (1996) 2693. 

[38] K.C. Jacobs, Spatially homogenous and Euclidean cosmological models with shear, 
Astrophys. J. 153 (1968) 661. 

[39] S. Hervik, Discrete Symmetries in Translation Invariant Cosmological Models, Gen. Rel. Grav. 33 
(2001) 2027 [arXiv:gr-qc/0105006]. 



- 22 - 



[40] C.B. Collins, Global structure of Kantowski-Sachs cosmological models, J.Math.Phys. 18 (1977) 2116. 

[41] J. Wainwright and G.F.R. Ellis, Dynamical systems in cosmology, Cambridge University Press, 
Cambridge - UK (1997). 

[42] C.B. Collins and J.M. Stewart, Qualitative cosmology, MNRAS 153 (1971) 419. 

[43] S. Hervik, W. C. Lim, P. Sandin and C. Uggla, Future asymptotics of tilted Bianchi type II 
cosmologies, Class.Quant.Grav. 27 (2010) 185006 [arXiv:gr-qc/1004.3661]. 

[44] A. Coley and S. Hervik, Bianchi models with vorticity: The type III bifurcation, Class. Quantum 
Grav. 25 (2008) 198001 [arXiv:gr-qc/0802.3629]. 

[45] L. Campanelli, P. Cea, G.L. Fogli and A. Marrone, Testing the isotropy of the universe with type la 
supernovae, Phys.Rev.D 83 (2011) 103503 [arXiv:astro-ph/1012.5596]. 

[46] J.D. Barrow, R. Juszkiewicz and D. H. Sonoda, Universal rotation: how large can it be?, 
MNRAS 213 (1985) 917. 

[47] T. Koivisto and D. F. Mota, Accelerating Cosmologies with an Anisotropic Equation of State, 
Astrophys.J. 679 (2008) 1 [arXiv:astro-ph/0707.0279]. 

[48] T. Koivisto and D. F. Mota, Anisotropic Dark Energy: Dynamics of Background and Perturbations, 
JCAP 0806 (2008) 018 [arXiv:astro-ph/0801.3676]. 

[49] T. S. Koivisto, D. F. Mota and C. Pitrou, Inflation from N-Forms and its stability, JHEP 0909 
(2009) 092 [arXiv:astro-ph/0903.4158]. 

[50] D. H. Lyth and A. R. Liddle, The Primordial Density Perturbation: Cosmology, Inflation and the 
Origin of Structure, Cambridge University Press, Cambridge - UK (2009). 



- 23 - 



Initial conditions 



D < 4 at time (a) 



Spacetime X Y Z Q, K V Q = 50 Q = 1000 



BI 


0.4 


-0.6 


0.5 






6.0 


5.7 


BI 


—0.4 


—1.1 


1.5 






13.4 


6.0 


BI 


—0.4 


0.3 


-0.7 






6.1 


5.8 


BI 


0.1 


0.7 


-0.7 






7.4 


5.8 


BI 


0.3 


-1.4 


-0.1 


- 


- 


9.5 


5.8 


BII 


-0.2 


-0.3 


0.5 


0.1 


-0.1 


2.9 


4.7 


BII 


-0.3 


-1.1 


-1.2 


0.2 


-0.4 


3.7 


5.2 


BII 


-0.1 


0.1 


-0.9 


0.1 


-0.1 


3.3 


4.8 


BII 


-0.5 


0.2 


-1.1 


0.1 


-0.1 


3.5 


5.2 


BII 


-0.6 


0.8 


-0.5 


0.3 


-0.2 


3.8 


5.0 


Bill 


-0.5 


-0.8 


0.6 


0.1 




8.1 


5.8 


Bill 


0.1 


0.8 


-0.4 


0.2 




7.4 


5.7 


Bill 


-0.1 


-1.5 


-1.3 


0.2 




12.7 


6.5 


Bill 


-0.6 


-1.7 


-0.3 


0.1 




14.1 


6.6 


Bin 


0.3 


-0.9 


0.5 


0.3 




8.2 


5.9 



Table 5. Simulations for Bianchi type spacetimes with initial conditions generated by a random 
number generator. D is the distance to fix-point (a) defined in (5.9). When D < ^ the universe is 
in the linear regime close to fix-point (a), and the distance to (a) is much smaller than the distance 
to (b). Two runs with different values for Q, while A = 0.1 in all runs. 



# 




Initial conditions 




D < — 


at time (a) 


H = at time (a) 


X 


Y 


Z 




Q = 50 


Q = 1000 


Q = 50 


Q = 1000 


1 


0.6 


-1.5 


2.2 


-0.8 


12.7 


6.4 






2 


0.8 


1.5 


-0.4 


-0.7 


12.4 


5.8 






3 


0.6 


-1.4 


0.3 


-0.5 


9.4 


0.0 






4 


-0.1 


2.1 


2.4 


-0.9 






0.3 


0.3 


5 


0.7 


-0.5 


2.0 


-0.9 


12.7 


6.1 






6 


2.2 


4.0 


3.9 


-9.3 






0.03 


0.03 


7 


1.6 


-2.5 


-0.2 


-3.1 






0.13 


0.13 


8 


1.5 


-4.9 


-3.8 


-8.7 






0.03 


0.04 


9 


-0.8 


-1.8 


2.6 


-2.3 






0.17 


0.17 


10 


-0.9 


2.7 


-4.9 


-5.1 






0.04 


0.04 


11 


0.6 


-1.4 


0.7 


-23 


3.8 


5.5 






12 


0.2 


-0.6 


-1.7 


-22 


4.1 


5.6 






13 


0.7 


2.0 


1.9 


-28 


5.2 


5.9 






14 


0.9 


1.9 


-0.1 


-21 


5.1 


5.5 






15 


-0.7 


-2.2 


-1.0 


-24 


4.5 


6.0 







Table 6. Simulations for KS with various initial conditions. D is the distance to fix-point (a) defined 
in (5.9). Two runs with different values for Q, while A = 0.1 in all runs. Initial values for where 
picked uniformly by a random number generator from the intervals (—1.0, —0.1), (—10.0, —2.0) 
and (—30, —20), respectively, for simulations #1-5, #5-10 and #11-15. In simulations #6-10 the 
potential energy of the scalar field is initially sub-dominant, while in #11-15 fiy ~ 
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Name Spacetime X Qk Stability 

(PI) FLRW 1 saddle 

(P2) BII |(l + 3w) ^(l + 2w-3u; 2 ) attractor 

(P3) Bill || attractor 

(P4) BI ±10 unstable 

Table 7. Post-inflationary fix-points. The unstable fix-point (P4) is either a saddle or a repeller 
depending on the spacetime and the sign of X (see eigenvalues). 



Name Spacetime Eigenvalues 

(PI) FLRW (l + 3w, -|(l-w)) 

(P2) BII (-|(1 -u), -|(1 -w)) 

(P3) Bill (-|, -3w) 

(P4) BI (4±2si, 3(1 -u)) 

Table 8. Real part of eigenvalues for post-inflationary fix-points. 
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